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Abstract 

We present a formulation of a matrix model which manifestly possesses the general co- 
ordinate invariance when we identify the large N matrices with differential operators. In 
order to build a matrix model which has the local Lorentz invariance, we investigate how 
the so(9, 1) Lorentz symmetry and the u{N) gauge symmetry are mixed together. We 
first analyze the bosonic part of the model, and we find that the Einstein gravity is re- 
produced in the classical low-energy limit. We then present a proposal to build a matrix 
model which has M = 2 SUSY and reduces to the type IIB supergravity in the classical 
low-energy limit. 



1 e-mail address : azuma@gauge.scphys.kyoto-u.ac.jp 

2 e-mail address : hkawai@gauge.scphys.kyoto-u.ac.jp 



1 Introduction 



A large N reduced model has been proposed as a nonperturbative formulation of superstring 
theory [00. It is defined by the following action: 

S = ~j2 Tr (\[A a ,A b }[A a ,A b } - ^T a [A a ,^ . (1.1) 



It is a large N reduced model ]T2|| of 10-dimensional M = 1 super Yang-Mills theory. Here, -0 



is a 10-dimensional Majorana-Weyl spinor field, and A a and ip are N x N hermitian matrices. 
This model is related to the type IIB superstring theory in that it is at the same time obtained 
by the matrix regularization of the Green-Schwarz action of the type IIB superstring theory 
(IIB matrix model is extensively reviewed in ||). 

IIB matrix model has several evidences of describing gravitational interaction; such as the 
exchange of gravitons and dilatons which appears in the one-loop effective Lagrangian[[|. How- 
ever, it is an interesting issue to pursue a matrix model which describes the gravity more 
manifestly. Our goal is to build a matrix model with manifest general coordinate invariance, 
which is the fundamental principle of general relativity. In order to achieve this goal, we should 
consider a matrix model with local Lorentz invariance. This issue has been investigated in the 



context of supermatrix model in [10| [IT 



In these papers, the matrix model based on the super Lie algebra u(l|16, 16)[|7|||[] has 
been investigated as a generalization of IIB matrix model to analyze how the local Lorentz 
invariance is realized in a matrix model. This model has the following two features. One is 
the extended Lorentz symmetry, in which the infinitesimal parameters of so(9, 1) are extended 
to the elements of u(N). In IIB matrix model, the information of the spacetime is embedded 
in the eigenvalues of the large N matrices, and the general coordinate invariance is related to 
the u(N) gauge symmetry Therefore, in order to realize the coordinate-dependent Lorentz 
symmetry, it is natural to regard the parameters of so(9, 1) as the elements of u(N). The other 
feature of the it(l|16, 16) supermatrix model is the existence of a higher rank tensor field. In 
order to formulate a matrix model in the curved spacetime, a spin connection term containing 
rank-3 gamma matrices must be included in the fermionic term. The it(l|16, 16) supermatrix 
model has been thus investigated as the model with higher rank tensor fields coupled to the 
fermions. 

In this paper, we inherit this idea, and attempt another formulation of a matrix model 
which manifestly possesses the local Lorentz invariance. The starting point is to identify the 



large N matrices with differential operators! 17]. This idea stems from the interpretation of the 



spacetime embedded in large N matrices. In the twisted reduced model (T^][|T4|, the matrix A a 
represents the covariant derivative in the spacetime, whereas in IIB matrix model the matrix A a 
itself represents the coordinate of the spacetime. The relationship between these two viewpoints 
is elucidated by expanding the matrix A a around the classical solution p a as 

A a =p a + a a , (1.2) 

where [p a ,Pb] = iBab and B a b is a real c-number. In |J, it has been pointed out that IIB 
matrix model is identified with the noncommutative Yang-Mills theory by this expansion. The 

3 so(9, 1) is a subalgebra of w(l|16, 16). For the detailed relationship between the super Lie algebra and the 
Lorentz symmetry, see (8| . 
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fermionic term of IIB matrix model i^TripT a [A a , ip] reduces to the action of the fermion in the 
flat space 



d°xi)(x)iT a {d a il){x) + [a a {x),^{x)]), (1.3) 

in the classical low-energy limit. 

With this idea in mind, we attempt to formulate a matrix model with the local Lorentz 
invariance, which in the classical low-energy limit reduces to the fermionic action on the curved 
spacetime 

S F = J d d xe(x)^(x)iT a e a i (x) [d^{x) + [A{x)^{x)\ + ^T bc u lbc (x)tP(x)^j . (1.4) 

Throughout this paper, the indices a, b, c, • • • and i, j, k, ■ • • both run over 0, 1, • • • , d — 1. The 
former and the latter denote the indices of the Minkowskian and the curved spacetime, respec- 
tivelyQ d is the spacetime dimension, and we focus on the case in which d = 10. In considering 
the correspondence with the matrix model, we need to absorb e(x) into the definition of the 
fermionic field, since it is not / d d xe(x) but / d d x that corresponds to the trace of the large iV 
matrices. And we regard the "spinor root density" 

vj/(x) = e ^(x)ip(x) (1.5) 

as the fundamental quantity. Then, the action is rewritten as 

S F = J d d x^(x)e^(x)iT a e a i (x) {di(e-^(x)^(x)) + [Ai(x),e-^(x)^(x)} 

+ ^T bc u ibc {x)e-^{x)m(x) 



d d x \^(x)iY a 



e a \x)di + ^e c \x)to lca (x) + e a \x)e^ (x)(^e a( x )) 



+^(x)r a eo i (x)[A i (x),*(x)] + ^(x)r aia2a3 e [ai i (x)a; ia2a 3 ] (x)*(x) 

d d x ^(x)iT a e a l (x)(d^(x) + [Ai(x),V(x)]) + l -^(x)T a ^ as e [ai l (x)u ia2Wi] (x)*(x) 

(1.6) 

where we have utilized in the last equality the fact that the fermionic field ty(x) is a Majorana 
one. The corresponding matrix model is formulated as 

S' F = ^Tr^T a [A a ,ij} + % -Tr^Y abc {A abc ^}. (1.7) 



4 In this paper, we adopt the following notation of the gamma matrices: 

{T a , r b } = 27f fc , where rf h = diag(-l, + !,•■•, +1). 



The gamma matrices are real and satisfy 



(r a y = ( T r a ) = | 



-r a (a = o) 

+r a (a = 1,2, ••■,9) 
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In promoting the action ( |1.6| ) to the matrix model, we have identified the covariant derivative 
with the commutator with the rank-1 matrix A a . The rank-3 term is a naive product of the spin 
connection and the fermion, and it is natural to promote the product to the anticommutator 
of the large N hermitian matrices^, i is necessary in the rank-3 term so that the action should 
be hermitian. The action (|1.7|) is an addition of the rank-3 term to the fermionic term of the 
action of IIB matrix model. Especially when if; is a Majorana fermion, ( [T77 ) is equivalent to 
the following action: 

S£ = Tnf>{T a A a + iT a ^A aia2a3 )ij. (1.8) 

Let us next consider the local Lorentz transformation of this matrix model. Originally the local 
Lorentz transformation of the fermionic field 8ip{x) = ^T aia2 s aia2 (x)if)(x) should be promoted 
to the anticommutator of the hermitian matrices as 5if> = ^T aia2 {e aia2 ,ip}. However, since 
it is difficult to find an action invariant under this transformation, we instead promote this 
transformation to the matrix version as 

Slf , = ^T a ^e aia2 ifj, (1.9) 

and we take the action ( |1.8j ). The hermiticity of the fermionic field is now sacrificed, and the 
actions (|1.7| ) and (|1.8| ) are no longer equivalent. The price we must pay for this alteration 
is that the product A a if) does not directly correspond to the covariant derivative. The local 
Lorentz transformation of the action (11.81) is 



SS'£ = ~^Tri> [T a A a + iT a ^A aia2a3J T b ^s blb2 ] if;. (1.10) 

Since we are now considering the local Lorentz transformation and their space-time dependent 
parameters are promoted to u(N) matrices, the infinitesimal parameters e a b are u(N) matrices. 
Then, the commutator 

r„"paia 2 a 3 a v bib2 F^ ,1 — 1 rr aia2a3 r bib2 M a p, , \ -u _ rr aia2a3 r bib2 \\ A , i 

I „ ' Z w ' 

rank-3 rank-1, 5 

(1.11) 

must include the rank-5 gamma matrices in 10 dimensions, and the action ( |1.8| ) is not invariant 
under the local Lorentz transformation. Repeating this procedure, we find that it is necessary 
to introduce all the terms of odd-rank gamma matrices to make the action invariant. 

The similar thing holds true of the generator of the local Lorentz transformation ( |1.9|) . In 
order for the algebra to close, only the rank-2 terms are not sufficient, because the commutator 

[J- Eaia 2 ,L £ bife 2 J~9 ] > J_ [ i Sa i a 2 ) £ b 1 b 2 J + 9 A \ I ^ ai ° 2 ' ^i^J \ l - LZ ) 

rank-2 rank-0, 4 



5 For the readers' convenience, we summarize the commutation relations of the hermitian and anti-hermitian 
operators. Let H and A be the set of the hermitian and anti-hermitian operators, respectively. When h, hi, h% € 
H and a. ai,ai 6 A, their commutation relations are as follows: 

[fti,/i 2 ] £A, [h,a]eH, [fli,fl 2 ]eA, {hi,h 2 }eH, {h,a}eA, {ai,a 2 }eH. 
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includes the rank-4 gamma matrices. Likewise, we find that the algebra of the local Lorentz 
transformation must include all the even-rank gamma matrices. 

Our formulation of the matrix model is based on this observation, and we analyze a matrix 
model which includes all odd ranks of 10-dimensional gamma matrices coupled to the fermion 
from the beginning. We identify these matrices with differential operators, and investigate 
whether our matrix model includes the supergravity in the classical low-energy limit. 

We mention the relation between our new proposal and the original IIB matrix model. In 
the quantum field theory, some different models which share the same symmetry are equivalent 
in the continuum limit, known as universality. We expect that the similar mechanism may hold 
true of large N matrix models and hence that various matrix models may have the same large 
N limit. Thus, there is a possibility that our new model is equivalent to IIB matrix model. We 
can expect that both IIB matrix model and our new model are equally authentic constructive 
definition of superstring theory in the large N limit. Even if IIB matrix model is not an eligible 
framework, we expect that our model may be closer to the true theory than IIB matrix model, 
because our model has larger symmetry. 

This paper is organized as follows. In Section 2, we analyze the bosonic term of the matrix 
model as the preliminary step, and we clarify that the bosonic part reduces to the Einstein 
gravity in the classical low-energy limit. In Section 3, we consider the model coupled to a 
fermionic field. We present a proposal to formulate a matrix model which has M = 2 SUSY 
algebra and reduces to the type IIB supergravity. Section 4 is devoted to the discussion and 
the outlook. 

2 Bosonic part of the matrix model 

Before going to the analysis of the supersymmetric matrix model, we first investigate the bosonic 
part: 



S = Tr 



NxN 



tr 32X 32V(m 2 ) . (2.1) 



The uppercase Tr and the lowercase tr denote the trace for the NxN and 32 x 32 matrices, 
respectively. The large iV matrix m consists of all odd-rank gamma matrices in 10 dimensions: 



7 1 ?' 1 

lit llt'd'- \ 2j '"'01101213 j-j"^!-^ ,/<, ai ... a7 l. g| "(-ai ■■■09 L ) 



(2.2) 



where m ai ... a2n _ 1 are large N hermitian matrices, m satisfies 

rV+r = m, (2.3) 

and thus the action ( |2.1| ) is hermitian. Odd powers of m must not appear in the action, since 
they would transform the fermion to that of the opposite chirality. The function V(m 2 ) is 
discussed in Section EO. 
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2.1 Identification of large N matrices with differential operators 

We identify the space of large iV matrices with that of the differential operators JT7|. By this 
identification, we can describe the differential operators on an arbitrary spin bundle over an 
arbitrary manifold in the continuum limit simultaneously, because they are embedded in the 
space of large iV matrices (some simple examples are provided in Appendix 0). In the following, 
from the big space of large N matrices, we pick up a subspace consisting of the differential 
operators over one manifold. We regard m as the differential operators over this manifold. 
Then, we analyze the effective theory of the fields appearing in the expansion of the differential 
operators (the explicit form is given later in (|2.6|) ). 

The space of the differential operator is infinite-dimensional, and in general the trace Tr for 
this space is divergent. But, as we elucidate later, we choose a function V(m 2 ) which decreases 
exponentially, and thus the trace is finite. We want to identify the dimensionless matrix m 
with something similar to the Dirac operator, and we need to introduce D as an extension of 
the Dirac operator in the curved space acting on the "spinor root density" Q1.5Q . D has the 
dimension [(length) -1 ], and m is expressed as 



i 



m = t 2 D, where (2-4) 
i 1 i 1 

n A F a _i A "pai(i2a3 a poi—as A pava7 j A -pa\---ag 

^ -^a 1 - T ^ ^i aia2a3 L ^j^a 0l ... a6 l ^|^i ai ... a7 ± -p gj ^ai-'-ag 1 > 

(2.5) 

where r has the dimension [(length) 2 ], r is similar to the Regge slope a' in string theory, and 
is introduced as a reference scale. This parameter r is not a cut-off parameter. V(m 2 ) is an 
exponentially decreasing function, and the damping factor is supplied by the action itself. Since 
V(m 2 ) is a function of the dimensionless quantity m, r represents a damping scale. When we 
approximate the differential operators by finite iV matrices, an ^-dependent ultraviolet cut-off 
naturally appears. When we take N to infinity, this cut-off becomes infinitely small. But the 
scale r is completely independent of this ultraviolet cut-off, and takes a constant value even in 
the large N limit. Thus, we can fairly take iV to infinity and identify the large N matrices with 
the differential operators, at least when we investigate the effective theory at tree level. In this 
sense, our model differs from the induced gravity. 

Aji---a 2 „_i i n = 1, 2, - • • , 5) are hermitian differential operators and are expanded by the 
number of the derivatives. Because of the hermiticity, they are expanded as the anticommutator 
of differential operators and real functions: 

oo jk 

^ax-oan-a = aax-a 2 „-i(z) + E ' ' ' <\ > ^'"^ax-a^x (x)}. (2.6) 

k=l Z 

The indices in the parentheses are symmetric, while the other indices are antisymmetric. Since 
D is an extension of the Dirac operator in the curved space, we find it natural to identify the 
function aj^(x) with the vielbein of the background metric. Now, the operator D is expanded 
as 



D = ei(x) 



i<„'(.r)r [di + h: bc u ibc {x) 



(x) 



+ (higher-rank terms) + (higher-derivative terms). (2-7) 
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The other coefficients a^ 1 '"' lk \ 1 ... a2n _ 1 (x) are regarded as the matter fields, and a simple dimen- 
sional analysis reveals that they have the dimension 

a (il - ifc) ai ...a 2 „_ 1 (x)~[(length)- 1+fc ]. (2.8) 



2.2 Local Lorentz invariance 

We next investigate the local Lorentz symmetry of this action. The local Lorentz transformation 
is 

5m = [m,e], where 

1 i 1 

e- _ ,>.. J -pa.ia2 j pOi-04 -pcif -ae 

c toy ~ 2| c iia2 c ai---<24 gj 1 c ai---a6 

Z 1 
pai— a8c- i -pai - aio /o q\ 

The large N matrix e satisfies 

rVr° = £ , (2.io) 

and represents the local Lorentz transformation, with e ai ... a2n (n = 0, • • ■ , 5) being hermitian. 
As we have observed in Section 1, all the even-rank gamma matrices are necessary in order for 
the algebra to close. The invariance of the action is verified as 



5S = Tr 



tr 



(2V'(m 2 )m[m,e])] = 0, (2.11) 



where V'(x) denotes V'(x) = 9V d ^ ■ The cyclic property of the trace still holds true of the space 
of the differential operators, if we assume that all the coefficients damp rapidly at infinity; i.e. 

lim a^-^ ai ... a2n _ 1 (x) = 0. (2.12) 

\x\^+oo 

Under this assumption, the trace of the commutator between the function and the differential 
operator is 

d d x{x\{d j a^\ 1 ... a2n _ 1 {x))\x) = f d d x(d j a^\ 1 ... a2n _ 1 (x))(x\x). (2.13) 



The surface term of this integral vanishes due to the assumption ( 2.12|) , and the cyclic property 
of the trace is ensured even for the trace of differential operators. 



2.3 Determination of V{m 2 ) 

We next discuss the condition V{m 2 ) should satisfy. In this section, we set the background 
metric to be flat for brevity. We require that the differential operator in the flat space ttlq = 
iT^T a d a should be a classical solution of this matrix model. To investigate this condition, we 
consider the Laplace transformation of the function V(u) 

rod 

V(u) = / dsg{s)e- su . (2.14) 
Jo 
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The Laplace transformation ( j2.14 ) enables us to relate the trace Tr[trV(m 2 )] with the Seeley 



de Witt expansion of the trace of e rD (T|] [|T^], which is the expansion around Wq = — r<9 a <9 a . 



The bosonic part of the action is thus expressed as 

Tr[trV(m 2 )} = / dsg{s)Tr[tre- srD } 
Jo 







dsg(s)[ I JJL[... + (st) i 2 A- 2 (x) + (sr) % U_i(ar) 




d d x 




(27rr)f L=_ oo V-'0 



d^(s) S -5 +fc ) rM fc (x) , (2.15) 



where runs over all the integers^. Ak(x) are the coefficients of the Seeley de Witt expansion of 
the trace Tr[tr(e~ srD2 )]. And there appear the coefficients At{x) for negative k, too, because 
the fields a^ 1 "" lk \ 1 ... a2n _ 1 (x) possess the dimension (|2.8| ). 

In order for m to be a classical solution, the linear term of the fluctuation around mo 
should vanish in the action. Since only the scalar fields can constitute a Lorentz invariant 
linear term, we focus on the linear terms with the indices contracted as a a ^ aiin "' ilil \x) (/ = 
0, 1,2, •••). Since a a ^'(x) is identified with the vielbein, the cancellation of the linear term 
Q>a { x ) means the absence of the cosmological constant in this matrix model. The linear terms 
of their derivatives (dj 1 ■ ■ ■ dj m a a ^ awi '" HH ^ {x)) disappear after integrating in the action. A simple 
dimensional analysis indicates that the linear terms of the fields 

OB («iii-w)( x ) are included in 

A-i{x). Therefore, we demand that the Seeley de Witt coefficients A-i(x) (/ = 0, 1,2, •••) 
should be canceled in the action. 

The condition for 7% to be a classical solution is thus translated into the following statement: 



o 



dsg(s)s-— n = 0, for n = 0, -1, -2, • • • . (2.16) 
Now, let us rewrite the condition ( 2.16| ) in terms of the function V(u). Noting 



poo foo roc poo 

/ duViuju"- 1 = du dsg{s)e- su u a - 1 = T(a) / dsg(s) S - a , (2.17) 
Jo Jo Jo Jo 

we discern that the function V{u) must satisfyf] 

duu* +n V(u) = 0, for n = -1, 0, 1, 2, • • • . (2.18) 



6 A dimensional ana lysis shows that the terms of the order 0(t for half integer I do not appear in this 

expansion. From (2.8), the dimension of the term nr=i(^ji '"®3p a k,i ^bi---b 2m -ii x )) ^ s [Q- en gth)' D ], where 



V = — n + X)™=i( — Pi + This term is included in the Seeley de Witt coefficient Av{x). On the other hand, 
the number of the indices in this term is S"=i(2wi — 1 + pi + Zj), which is equal to T> modulo 2. When V 
is odd, the total number of the indices is also odd and this term cannot contract the indices to constitute a 
scalar. Therefore, the coefficient Ai(x) for half integer I is prohibited due to the Lorentz invariance of the action. 
However, they may survive when d is odd, with the antisymmetric tensor e^...^ being accompanied. 

7 Sincc | + n > for n = —1, 0, 1, 2, • • •, we need not worry about the possibility that T(4 + n + l) vanishes. 
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There are various choices for such functions, and one example is the following function:^ 

gi- 1 (exp(-^i)sin(Mi)) 
Vo{u) = —j— . (2.19) 

Going back to the general case, when V(m 2 ) is chosen in order that m = ir^T a d a becomes a 
classical solution, we understand that the action reduces to the Einstein gravity in the classical 
low-energy limit. The absence of the linear term 0( a ) a (x) ensures that the graviton is massless, 
since the mass terms (including the cross-terms with the other fields) vanish due to the general 
coordinate invariance. Then, the Einstein gravity is derived from the Seeley de Witt expansion 
of the trace Tr[tr(e~ rD )] when we bequeath only the curved space Dirac operator in the 
expansion fl2.7|) fl 

Tr[tr(e- TL>2 )} = [ d d x—^e(x) (r^- + ■■), (2.20) 



(2ttt) 



2 



where 32 comes from the trace £rl32 X 32- 

It is obvious from the dimensional analysis that this term is included in A\(x). When A±(x) 
survives in the action, i.e. / °° dsg(s)s~2 +1 ^ or equivalently J °° duV{u)u^~ 2 ^ 0, the trace 
Tr[trV(m 2 )] includes the Einstein gravity in the Seeley de Witt expansion. 

We next investigate the mass terms and the kinetic terms of the matter fields. A dimensional 
analysis shows that the mass terms and the kinetic terms (including the cross terms) are 
included in the following Seeley de Witt coefficients: 

mass terms: a ( * 1 "'* fe) ai ... a2n _ 1 (x)a (jl "^ ) ai ... 02 „_ 1 (x) G A x _k+i{x), (2.21) 
kinetic terms: (^ 1 a (il "- i * ) ar .. a3B _ 1 (a;))(^ a a (, ' 1 - J ' ) ai ... aa „_ 1 (a;)) G A 2 _k±i(x), (2.22) 

where k + I is an even number. In particular, the mass terms and the kinetic terms of the odd- 
rank anti-symmetric fields a 0l ... a2n _ 1 (x) are included in Ai(x) and A2(x), respectively. Generi- 
cally0, these fields are massive and are decoupled in the classical low-energy limit. 

Next, we consider the fields ai-aa n -i( x )- As we w * n explain in Section 3, especially the 
fields a^'iai-a,2n{ x ) ( w ith n = 1,2, 3, 4) are identified with the anti-symmetric tensor fields in the 
type IIB supergravity. Their mass terms are included in Ao(x), whereas the kinetic terms reside 
in Ai(x). These fields are therefore massless. This observation provides a positive evidence in 
constructing a matrix model which reduces to the type IIB supergravity, which is the goal of 
the discussion of the next section. 



We can verify that this function satisfies ( 2.18| ) as follows [fl6| . We note that J °° dyy rn e ay 



for a = cxp(^) = and m = 0, 1, 2, • • -. This is a real number when m — 3 is a multiple of 4. Taking the 
imaginary part of the both hand sides, we obtain dyy 4n+3 sin(^=) exp(— -2*=) = 0, for n = 0, 1,2, We 



make a substitution u = ^- and integrate by part to verify that ( 2.19 ) satisfies ( 2.18 ). Another choice is the 



function V(u) = — — ^"y , which corresponds to the limit a = exp(^). In this case, we need to introduce a 
proper damping factor. 

9 The square density e 5 (x) can be treated by using the cyclic rule of the trace. 

10 The computation of the trace Tr\tr(eT TD )] indicates that the mass term a a (x)a a (x) vanishes, but a a (x) 
is presumably massive because there is no reason that the cross-terms a a (x)a^ 1 '"^ 2 '^ a {x) also vanish. 
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However, it is not clear whether the higher-spin fields a ( - n "" lk \ 1 ... a2n _ 1 (x) (k = 2, 3, • • •) are 
massive, since the mass terms and the kinetic terms belong to the coefficients A\-k(x) and 
A2-k{x) respectively and both of them vanish in the action. 

3 Supersymmetric Action 

In the previous section, we have investigated the bosonic part of the matrix model, and we 
have elucidated that the Einstein gravity is derived from the bosonic action in the classical 
low-energy limit. 

We next add the fermionic term and present a proposal to construct a matrix model which 
has M = 2 SUSY and reduces to the type IIB supergravity in the classical low-energy limit. If 
we succeed in building such a matrix model, this will be an extension of IIB matrix model with 
the gravity encoded more manifestly and provide a strong evidence that a matrix model is an 
eligible framework to describe the gravitational interaction. The suggestion in this section is 
rather conjectural. We start with the following action: 

Ss = Tr[trVs(m 2 )] + Tripmip. (3.1) 

Here, m is a large N matrix already defined in (|2.2|) . m is identified with the differential 
operators and expanded by the number of the derivatives in the same way as is explained in 
the previous section. 

ijj is a Weyl fermion, but is no longer hermitian. As we have observed in Section 1, it is 
difficult to construct a matrix model which is invariant under the local Lorentz transformation 
of the type Sip = ^T ab {e a b, ip}. Then, we define the local Lorentz transformation of the matrix 
model as 5ip = sip and sacrifice the hermiticity of the fermionic field. This fermionic field is 
also identified with differential operators and expanded as 

^(xW + Ex^W^-^e-^, (3.2) 

where the fermionic fields x^ l ' H \ x ) possess the dimension [(length)']. We temporarily intro- 
duce the damping factor e~^ rD so that the trace Trtprmp should be finite. The power a is 
chosen in accordance with the function Vs(m 2 ) and, for example, when we choose the function 
(|2.19[ ), the power is a = |. 

This action is invariant under the local Lorentz transformation 

Sm=[m,e], Sip = eip, (3.3) 

where e is already defined in (|2.9|) . The transformation of ip is readily seen to be Sip = —ipe. 
The local Lorentz invariance is verified as 

SS S = 2Tr[tr(Vs(m 2 )m[m, e))) + Tr[tr(i>[m, e]ifi)] = 0, (3.4) 

when we assume that the coefficients Or il '" 1,k ai—aa n -i( x ) anc ^ X^ ll ' lk \ x ) damp rapidly at infinity. 

The function Vs(m 2 ) is to be determined so that it satisfies the following conditions. First, 
the differential operator in the flat space m = iT^T a d a must be a classical solution, and 
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Vsijn 2 ) must satisfy the criterion fl2.18| ). Second, in order for the model to reduce to the type 
IIB supergravity in the classical low-energy limit, only the fields aS l \ aia2 ... a2n (x) (n = 1, 2, 3, 4), 
which are interpreted as the anti-symmetric even-rank fields in the type IIB supergravity, must 
be massless and the other fields must be massive. It is not clear whether ( |2.19| ) is an eligible 
function, because we have yet to understand whether the higher-spin fields are massive. This 
problem will be reported in more detail elsewhere. 

We next consider the structure of the supersymmetry of this matrix model. The SUSY 
transformation is given by 

5 e m = eip + ipe, 5 e ip = 2V' s {m 2 )e. (3.5) 

It readily follows that the SUSY transformation of ip is 

5 e ip = 2eVg(m 2 ). (3.6) 

The invariance of the action under this SUSY transformation is verified as 

5 e S s = Tr [tr ({2V' s {m 2 )m(eil) + ipe)\ + ^(<a/J + ipe)ip 

+2ifjmVs(m 2 )e + 2emV^(m 2 )^)] =0. (3.7) 

We analyze the commutator of this SUSY transformation especially when it is possible to 
perform the Taylor expansion of Vs(u) around it = as 

oo 

Vs(u) = £ (3-8) 
fe=i ZK 

in order to simplify the analysis using the equations of motion 

= 2mi, = 0, ^ = 2^jm = 0. (3.9) 
dip dip 

The following analysis is not applied to the function (|2.19|) , because the Taylor expansion 
around u = is impossible^. However, in the case of (|3.8|), the commutator can be simplified 

as 

[<^]TO = 2[£e-e£,^(m 2 )], (3.10) 

R , W = 2, (t SHr*)- ™ _ ^)- Vm \ _ (3 . n) 
y m z m A J 

The explicit derivation is presented in Appendix. |B|. 

In order to see the structure of the M = 2 SUSY, we separate the SUSY parameters into 
the hermitian and the antihermitian parts as 

e = ei+ie 2 , £ = & + i£ 2 , (3.12) 



ii 



The invariance of the action under the local Lorentz transformation (^4) and the SUSY transformation 
( 3.7) holds true without assuming the Taylor expansion around u = 0. 
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where e±, e 2 , £1 and £2 are Majorana-Weyl fermions. For simplicity, we assume that e and £ are 
c-numbers and that the background metric is flat. We now clarify that the translation of A a 
comes from the quartic term of m in the Taylor expansion (|3.8|) : 

[S e ,Ss}A a = ±tr([5 e ,8s}mT a ) 

= — Y a 2k tr(iem 2k - 2 Y a - e^m 2k ' 2 T a - m 2fc ~ 2 £er a + m 2fc - 2 e£r a ) 
1 00 

= £ Me> 2fc ~ 2 , Tie - e[m 2k - 2 , T a ]0 

1,3 k=2 

= ^(£[1^, T a ]e - e[T b ^ b \T a }OA bl A b2 + ■■■ = ^(£Pe - eV^A, A a ] + ■ ■ ■ 
lb lb 

= jfaVei + ^re 2 )[A u A a } + ---, (3.13) 

where we have utilized the fact that the SUSY parameters e and £ are c-numbers and thus 
that the matrix m commutes with the SUSY parameters. Let us concentrate on the SUSY 
transformation of the field a a (x). The commutator [A{, A a ] represents the translation of a a (x) 
and its gauge transformation: 

[Ai, A a ] = [idi + a,i(x),id a + a a (x)) H 

= i{diaa{x))^ -i{d a ai{x)) + [a>i(x), a a {x)} | H • (3.14) 

translation gauge transformation 

Therefore, we discern that the commutator of the SUSY transformations constitutes a bona 
fide translation for the vector fields. 

However, this is not the case with the fermionic fields. The commutator is computed as 

n 

[S e , = -£ a 2k ^m 2k ~ 3 e - em 2k - 3 + ■■■ = -a 4 (£Fe - eT^A, + . . . 

k=2 

= -2a 4 (6r j e 1 + l 2 Y j e 2 )^A j + ■■■. (3.15) 
We explore the term ipAi more carefully: 

iPAj = #9, + • • • = ^(1)9, + f A (!l - !,) (^)3n ■ • • 9,9, j e"^ 2 )" + ■ • • , (3.16) 

where • • ■ denotes the omission of the non-linear terms of the fields. Therefore, each fermionic 
field is transformed as 

[S € ,S 6 ]x(x) = + ---, 

[8 e , <yx (<1 *"* ,+l) (z) = -204(^1^6! + &Pe 2 ) X {{il ~ il) (x)8 i ^ j + ■■■. (3.17) 



The commutator of the SUSY transformation constitutes a translation for the vector fields 



while the same does not hold true of the fermionic fields. However, the analysis (|3.10|) — ( 3.1 7|) 



is applicable only to the case in which we can perform a Taylor expansion of Vs(u) around 
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u = 0. We speculate that we may be able to construct a well-defined model with Af = 2 SUSY 
for the functions like ( |2.19|) , which do not have the Taylor expansion around u = 0. 

We surmise that the fermionic fields x( x ) an d X i x )> which are identified with the dilatino 
and the gravitino respectively, become massless due to the supersymmetry, when the even- 
rank fields oy\ai-a2 n ( x ) are massless. If we choose the function Vs(m 2 ) properly, the classical 
low-energy limit of this matrix model would reduce to the type IIB supergravity, and it is an 
interesting future problem to investigate this conjecture more deeply. 

4 Conclusion 

We have hitherto investigated the possibility to realize the general coordinate invariance mani- 
festly in a matrix model. And in order to build a matrix model which describes the gravity, we 
should consider the local Lorentz invariance of the model. In realizing the local Lorentz invari- 
ance in a matrix model, the parameters of the Lorentz symmetry must be promoted to u(N) 
matrices, because the information of the spacetime is encoded in the eigenvalues of the large N 
matrices. The model thus has to include all the odd-rank gamma matrices in 10 dimensions. 
The starting point to construct such a matrix model is to identify the large N matrices with 
differential operators. By this identification, it is possible to describe the differential operators 
on an arbitrary spin bundle over an arbitrary manifold simultaneously in the space of large N 
matrices. 

We have first considered the bosonic part of the model, and clarified that it reduces to the 
Einstein gravity in the classical low-energy limit. And we have extended the observation to 
the supersymmetric action, and have presented a proposal to construct a matrix model which 
has Af = 2 SUSY and reduces to the type IIB supergravity. It is a future problem to seek the 
explicit form of the function Vs(m 2 ). If we find a proper model, this would be a nonperturbative 
formulation of the type IIB superstring theory, and become an extension of IIB matrix model 
with the information of the gravity encoded more manifestly. It is also interesting to elucidate 
the relation between our model and IIB matrix model. More analysis of these issues will be 
reported elsewhere. 
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A Examples of the differential operators in the space of 
large N matrices 

This appendix is devoted to some simple examples of the differential operators embedded in the 
space of large N matrices. The space of the large N matrices includes the differential operators 
on an arbitrary spin bundle over an arbitrary manifold simultaneously. We now consider the 
differential operators of the scalar field on two different bundles over Si. 
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Figure 1: The way the spaces of the differential operators are embedded in the space of large 
N matrices. 



1_ 
N 



_2 
N 



N_ 
N 



(1) trivial bundle with the periodic 
condition f(1)=f(0). 



(2) Z2-twisted bundle with the antiperiodic 
condition f(1)=-f(0). 



Figure 2: Differential operators on (l)the trivial bundle and (2)the Z 2 -twisted bundle over £1 



(1) trivial bundle 

We first consider the trivial bundle with the periodic condition f(l) = f(0). We discritize 
the region < x < 1 into small slices of spacing e = Then, the differential operator is 
approximated by the finite difference as 



(T) . 1 /7(¥)-/(#) , /(#)-/(¥)' 



N 
~2 



k + 1 



f 



k-1 



2 \ e ' e J -1 V" V v / V - v 

Due to the periodic condition, this finite difference is expressed by the large iV matrix as 

-1 \ 

-1 1 



• (A.I) 



d x 



( 1 
-1 1 

-1 1 



V 



-1 



(A.2) 



/ 



(2)Z 2 -twisted bundle (Mobius strip) 

We consider the similar problem with respect to the Z 2 -twisted bundle, in which the antiperiodic 
condition f(l) = — /(0) is imposed. Paying attention to this antiperiodicity, we understand 
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that the finite difference in the discritized space is expressed as 

1 \ 



d x ^ A = — 
2 



/ 1 
-1 1 

-1 1 



(A.3) 



-1 0) 



B Proof of ( gTtg ) and 



This appendix is devoted to the explicit computation of the commutation relations of the SUSY 
transformation. First, ( |3.1U| ) is verified by considering the differences of the following two SUSY 
t r ansf or mat ions : 

m m + + $1 h m + (e + + -0(e + £) + 2£eV s (m 2 ) + 2Vg(m 2 )e£, 
m h m + eifi + ipe ^ m + (e + + -0(e + £) + 2e^V s (m 2 ) + 2V s {m 2 )£e. 



1> 



v 



( |3.11D is verified likewise, but we need to use the Taylor expansion of Vs(u) around u = 0: 

St 



Se 



5 e 



il> + 2V' 8 {m*)S 

oo 

i/> + 2V s (m 2 ){e + V) + E a 2k[(e^ + </>e 



rn 



2k-3 



m(e^ + 4>e)m 2k ~ 4 -\ h m 2k ~ 3 (etp + ipe)]£ 



k=2 



i) + 2V s {m 2 )e 



ip + 2V s {m 2 ){e + ip) + J2 a 2k[{& + ^)m 2k ~ 3 + + ^)m 2fc " 4 + • • • + m 2fc " 3 (^ + ^)]e. 



k=2 



Utilizing the equations of motion fl3.9|) and the fact that the combinations of the fermionic 
fields ipe and ip^ vanish because the fermions are Weyl ones, we find that the commutator of 
the SUSY transformation is 



,2k- 



[S e , S^ip = a 2k (ijem 2k 3 £ - 4>£m 

k=2 \ 



3 n nl ( Vst™ 2 ) - Vl(0) , - V s {m 2 ) - V^(0) 

3 e) = 2ip em— ; - sy ' £ - £m— — 1 - sy < 

\ m 2 m z 
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